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ABSTRACT 

The discovery of over 200 extrasolar planets with the radial velocity (RV) tech- 
nique has revealed that many giant planets have large eccentricities, in striking 
contrast with most of the planets in the solar system and prior theories of planet 
formation. The realization that many giant planets have large eccentricities raises 
a fundamental question: "Do terrestrial-size planets of other stars typically have 
significantly eccentric orbits or nearly circular orbits like the Earth?" Here, we 
demonstrate that photometric observations of transiting planets could be used to 
characterize the orbital eccentricities for individual transiting planets, as well the 
eccentricity distribution for various populations of transiting planets (e.g., those 
with a certain range of orbital periods or physical sizes). Such characterizations 
can provide valuable constraints on theories for the excitation of eccentricities 
and tidal dissipation. We outline the future prospects of the technique given the 
exciting prospects for future transit searches, such as those to be carried out by 
the CoRoT and Kepler missions. 
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1. Introduction 

Theorists liave proposed numerous mechanisms that could excite orbital eccentricities. 
Some of these mechanisms are expected to affect all planets independent of their mass (e.g., 
perturbations by binary companions, passing stars, or stellar jets; e.g., Holman et al. 1997; 
Laughlin & Adams 1998; Ford et al. 2000; Zakamska & Tremaine 2004; Namouni 2007), 
while the efficiency of other mechanisms would vary with planet mass (e.g., planet-disk or 
planet-planet interactions; e.g., Artymowicz 1992; Goldreich & Sari 2003; Chatterjee et al. 
2007). If the mechanism(s) exciting eccentricities of the known giant planets also affect 
terrestrial planets, then Earth-mass planets on nearly circular orbits could be quite rare. 
On the other hand, if large eccentricities are common only in systems with massive giant 
planets and/or very massive disks, then there may be an abundance of planetary systems 
with terrestrial planets on low eccentricity orbits (Beer et al. 2004). Thus, understanding the 
eccentricity distribution of terrestrial planets could provide empirical constraints for planet 
formation theories (e.g.. Ford & Rasio 2007) and shed light on the processes that determined 
the eccentricity evolution in our solar system. Since the discovery of transiting giant planets 
with eccentric orbits, authors have begun to consider the implications of eccentricities for 
transiting planets (e.g., Barnes 2007; Burke 2008). 

The CoRoT and Kepler space missions are expected to detect many transiting planets 
and measure their sizes and orbital periods, including some in or near the "habitable zone" 
(e.g., Kasting et al. 1993). The Kepler mission aims to determine the frequency of Earth- like 
planets, and study how the frequency and properties of planets correlates with the properties 
of their host stars (Basri et al. 2005). Since a significant eccentricity would cause the stellar 
flux incident on the planet's surface to vary, a planet's eccentricity affects its climate (i.e., 
equilibrium temperature, amplitude of seasonal variability) and potentially its habitability 
(Williams et al. 2002; Gaidos & Williams 2004). Our method could be applied to these 
planets to determine the frequency of terrestrial planets that could also have Earth-like 
chmates, and thus influence the design of future space missions that will attempt to detect 
and characterize nearly Earth-like planets (e.g.. Space Interferometry Mission-Planet Quest) 
and search them for signs of life (e.g.. Terrestrial Planet Finder). 

In ^ we show how the duration of a transit is affected by a planet's orbital eccentricity. 
We outline how to interpret the transit duration for transiting planets with both low signal- 
to- noise light curves ( §2.1.ip and high signal-to-noise ( §2.1.21) . We compare the magnitude of 
the effect on the transit duration to the expected precision of eccentricity constraints based on 
Kepler photometric data ( §2.2p and also the typical accuracy of stellar parameters ( §2.3p . We 
demonstrate that Kepler observations could be used to characterize the orbital eccentricities 
of terrestrial planets. We describe statistical approaches for analyzing the distribution of 
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transit durations of a population of transiting planets in §3] In §3.4[ we discuss the role 
of radial velocity observations for constraining eccentricities of Earth-like transiting planets. 
In §11 we conclude with a discussion of how the results could contribute to understanding 
the formation and evolution of terrestrial planets and address fundamental questions, such 
as "What is the frequency of terrestrial planets that have Earth-like eccentricities?" and 
"What is the frequency of terrestrial planets that pass through the habitable zone?" 



The total transit duration is defined as the time interval between the first and fourth 
points of contact (Fig. [1]). In the approximation that the mean planet-star separation (a) 
is much greater than both the stellar radius (i?^) and the planetary radius (Rp), the transit 
duration is much less than the orbital period (P), and thus the planet-star separation during 
the transit (dt) is nearly constant. Using these approximations, the total transit duration 
(tz)) is given by 



where e is the orbital eccentricity, r = Rp/ i?^ is the ratio of the planet radius to stellar radius, 
b = dt cos i/Ri, is the impact parameter, and i is the orbital inclination measured relative 
to the plane of the sky (Tingley & Sackett 2005). The planet-star separation at the time of 
transit is given by dt = a{l — ecosEt) = a(l — e^)/(l-|-ecosTj) = a(l — e^)/(l+ecosti;), where 
Tt and Et are the true and eccentric anomalies at the time of transit, and u is the argument of 
periastron measured relative to the line of sight. (Note that this differs from the convention 
for radial velocity determinations.) Thus, the ratio of the actual transit duration to the 
transit duration for the same planet on a circular orbit (assuming other transit parameters 
such as size and impact parameter are held fixed) is (1 + ecosuj)/y^l — e^). We have 
performed numerical integrations to verify that this approximation is accurate to better 
than ~ 0.1%, for eccentricities as high as 0.95 (Fig. [2]). 

We define the variable Tb = dt/{ayl — e^) to be the ratio of the transit duration of 
one planet to the transit duration for an identical planet with the same orbital period and 
impact parameter, but on a circular orbit. In general, transiting planets could have significant 
eccentricities, in which case Tb will deviate from unity (e.g., Tb ~ 1.8 for HAT-P2b; Bakos 
et al. 2007). For planets on low eccentricity orbits, \Tb — 1| ~ 0(e). For large eccentricities, 
Tb can range from zero to several (Fig. [2]). For the case of a central transit of an Earth-like 



2. Light Curve of an Eccentric Transiting Planet 
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planet orbiting a solar- mass star on a circular orbit, 

«„.13hou.(i^)'"(|.) [M^y\^,,r^'mr,rf-,Y\. (2) 

where is the stellar mass and fi = mp/M^ is the planet-star mass ratio. The planet-star 
mass ratio (/x = rUp/M^^) is typically negligible for planetary companions. (If it is large, 
then it could be measured by obtaining RV observations.) By assuming that the pericenter 
direction is randomly oriented relative to the line of sight, we can compute the probability 
distribution for as a function of eccentricity (Fig. [3]). 

2.1. Analysis of Transiting Planet Light Curves 

Here we outline the basics of how photometric light curves can be used to constrain the 
eccentricities of transiting planets. We discuss two limiting regimes: 1) low signal-to-noise 
(S/N) light curves that provide no constraint on the impact parameter, and 2) high S/N 
light curves that measure the impact parameter. 

2.1.1. Impact Parameter Not Measured (Low S/N) 

For any transit detection, photometric observations alone directly measure the orbital 
period (P), the total transit duration (tz)), and the depth of the transit. For the purposes 
of providing analytic estimates, we neglect the effects of limb darkening, so that the transit 
depth is given by the planet-star area ratio (r^). For some faint stars and/or small planets, 
it may be difficult to measure additional light curve parameters such as the ingress duration 
and the impact parameter. For these cases, the total transit duration, to, can be compared 
to t/) o, the transit duration expected for the same planet with the same orbital period (P) 
and mean stellar density (p*), but assuming a circular orbit and central transit (6 = 0). Even 
for transiting planets with low S/N light curves, photometric observations can measure the 
ratio 

^^llL= ( \ I \f^+^^A ^ ( tj, \ /yr A V3 (p^\ (1+p)^/^ 

^° tD,o \aVl^)\ 1 + ^ ] VlShoursy* Vpy Vpo/ (1 + 

(3) 

Since = rrip/M^ <^ 1 for planetary mass companions, a value of Tq significantly greater 
than unity can only arise for an eccentric planet, but a value of less than unity could be 
due to either a non-central transit {b > 0) or a non-zero eccentricity (e > 0). Thus, it is not 
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possible to measure the eccentricity for an individual planet with only low S/N photometry. 
However, it is still possible to characterize the distribution of eccentricities for a population 
of planets based on the observed distribution of Tq (see §3.2p . Since the above analysis can be 
applied to relatively faint stars with low S/N transit light curves, we expect that Kepler will 
discover many planets that can be included in statistical analyses of Tq (Basri et al. 2005). We 
note the actual analysis should consider the effects of limb darkening, which cause the depth 
of a high-latitude transit relative to be greater than a than an equatorial transit even for 
the same star-planet area ratio. While limb darkening precludes simple analytic expressions, 
such effects can be incorporated into Monte Carlo simulations. 

In ^ we will show that Kepler observations could constrain the distribution of eccen- 
tricities and reject plausible eccentricity distributions for terrestrial planets. The precision of 
the eccentricity constraints will depend on the precision and accuracy of the measurements 
of To. In §2.21 we show that the the orbital period, transit depth, and transit duration can 
be measured with sufficient precision that the accuracy for the measured fg will typically be 
limited by the deviation of the inferred stellar density (p^) from the actual stellar density 
(p^,), i.e.. To ~ To{pi,/ pi,y^^. The mean stellar density could be estimated using spectroscopy 
and stellar modeling ( §2.3p . resulting in a typical accuracy for Tq of ~ 5 — 15% (e.g.. Ford et 
al. 1999; Fischer & Valenti 2005; Takeda et al. 2007). The spread of transit durations due to 
the unknown inclination will have a negligible effect on the ability of low S/N transit light 
curves to constrain the eccentricity distribution. 



2.1.2. Impact Parameter Measured (High S/N) 

For high-quality light curves, photometric observations can measure the time of each 
point of contact and determine both the total transit duration (t^^) and tp (Fig. [1]). If we 
neglect limb darkening, then the impact parameter can be determined by the photometric 
observables, 6^ = [(1 — r)^ — 7^ (1 + r)^] / (1 — 7^) , where we define 7^ = tj^/tj-, (Seager & 
Mallen-Ornelas 2003). 

In practice, the transit shape will be affected by limb darkening (Fig. H]) and we must 
fit for the limb darkening parameters to determine the impact parameter (and planet-star 
radius ratio). Previous experience analyzing transit photometry of giant planets has found 
that the limb darkening parameters can be accurately estimated with stellar models and the 
uncertainty in limb darkening parameters typically introduces relatively modest additional 
uncertainties and have only small correlations with the other model parameters of interest 
(e.g.. Brown et al. 2001; Holman et al. 2006; Burke et al. 2007; Knutson et al. 2007; Winn 
et al. 2007; Torres et al. 2008). In order to verify that the uncertainty in limb darkening 
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coefficients will introduce only modest additional uncertainties, we have generated simulated 
Kepler V-band light curves for terrestrial planets transiting a. V = 12 magnitude solar-like 
star, including a quadratic limb darkening model (Mandell & Agol 2002; Claret 2000). We 
then calculate Ax^, varying the transit time, transit duration, the planet size, impact pa- 
rameter, and quadratic limb darkening coefficients (see Fig. [5]). While these are some modest 
correlations between the transit duration and limb darkening parameters, these correlations 
are a much smaller effect than the uncertainty due to the correlation between the transit 
duration and the impact parameter. Therefore, we conclude that the uncertainty in impact 
parameter will be the limiting factor and we can neglect the effects of limb darkening when 
estimating the power of transit durations for constraining the eccentricities of transiting 
planets (see §3.2p . Once the host star and transit time are known from Kepler observations, 
a future transit can be observed using ground-based observatories at multiple wavelengths. 
By combining the inferred stellar properties with limb darkening profiles from stellar models, 
one can determine the impact parameter from even relatively low precision data (e.g., Jha 
et al. 2000). Infrared observations would be particularly useful, as they minimize the effects 
of limb darkening and aid in the measurement of the impact parameter. Since the impact 
parameter is the same at all wavelengths, ground-based transit light curves could effectively 
eliminate the need to fit for this parameter when analyzing the Kepler photometry. Here, 
we assume that limb-darkening parameters can be well constrained by some combination of 
stellar modeling and external observations. 

For systems with high S/N photometry, the total transit duration, t^i, can be compared 
to t£, h, the transit duration expected for the same planet and star with the same observed 
orbital period, impact parameter (6), and mean stellar density, but assuming a circular orbit. 
Then, Kepler can measure the ratio 

The fact that the planet is transiting the star places a constraint on a combination of true 
anomaly and the direction of pericenter. Therefore, the actual value of Tb depends only on 
the eccentricity and direction of pericenter, and can be measured using a combination of 
transit photometry and stellar parameters (e.g., from spectroscopy and stellar modeling). 
Using high S/N photometry, we can accurately measure the ffist three terms on the right 
hand side of Eqn. |H The last term containing the planet-star mass ratio (/i) is negligible for 
planetary companions. Therefore, the accuracy of the measured Tf, will typically be limited 
by the uncertainty in the cube root of the mean stellar density {a 1/3). We can detect a 

non-zero eccentricity for an individual planet when e > / pj — 0.05 — 0.15 ( §2.31) . This 
is quite significant given the mean eccentricity of planets discovered by RV surveys is ^ 0.3 
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(ignoring planets with P < 10 days that may have been influenced by tidal circularization; 
Butler et al. 2006). 

2.2. Expected Precision of Kepler Measurements 

Here we present estimates of the timing precision possible with Kepler observations, 
assuming a flux measurement precision of aph ~400ppm during a one minute integration on 
a Y—12 star and uncertainties that scale with the square root of the photon count (Basri et 
al. 2005). For the sake of deriving approximate analytic expressions, we assume uncorrelated 
Gaussian uncertainties and ignore complications due to limb darkening in this paper 

To determine the precision of eccentricity constraints, we must estimate the precision 
of the constraints from transit photometry. The fractional transit depth (r^ in the absense 
of limb darkening) can be measured with a precision, (7^2 ~ aph/ \/tpAN^, where A is the 
rate of photometric measurements with precision Cph and Ntr — 4(yr/P) is the number of 
transits observed during the four year mission lifetime. Note that in this approximation 
the precision does not depend on the integration time, provided that it is significantly less 
than the ingress/egress time. In practice, an increased rate of measurements is valuable 
for constraining any brightness variations across the stellar disk (e.g., limb darkening, star 
spots, plage). 

The time of each ingress/egress {tin/out) can be measured with a precision o't-,^/^^^ — 
Oph^J M^injoutl ^ {R*/ Rpf 1 where ^tin/cmt is the duration of ingress/egress (Ford & Gaudi 
2006). Therefore, the sidereal orbital period can be measured with a precision of up ~ 
^tin/outl '\/'^ (^tr ~ 1) s-nd the duration of the transit can be measured with a precision of 
(7t^ = crti„/out a/VA^. Once there are at least two transits, the uncertainty in will domi- 
nate. 

The duration of ingress (Atj„) or egress(Atout) is given by Atin/out = e^D, where 
e(l — e) = r/ [(1 + r)^ — 6^] . For non-grazing transits, e <S 1, so the ingress/egress du- 
ration is 

i_ 1 
AU„ Atout ^ 7.1minutes (^^^ (^^^ ' ((1 + rf - b')^ (^^^ ' (1 + fx)^ r,. (5) 

Inserting fiducial values for a target star with an apparent magnitude of V, we expect a 
fractional precision for the mean transit duration of an individual planet to be 

^ e)^ © s)' '-"^-'^ ('-^^^ 

(6) 
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For a Jupiter-sized planet with a similar orbit and host star, the expected fiducial fractional 
uncertainty for Kepler observations decreases to ^ 3 x 10"'*. Ground-based observatories 
have achieved a precision of crt^/t_D — 1-5% for short-period giant planets (e.g., Holman et 
al. 2007). 

For modest eccentricities, Kepler's expected measurement precision will place a lower 
limit on the uncertainty of a given planet's eccentricity (cTe) of order ~ dj^/ti). Even for 
Earth-sized planets around V=14 stars, Kepler is expected to achieve the photometric pre- 
cision necessary to measure eccentricities as low as ~ 0.03. 

For the sake of completeness, we also estimate one additional timescale. For an eccentric 
orbit, there is a slight change in df between ingress and egress that results in a difference 
between the ingress and egress durations, 

/p\-i/3 /^^x /p \ /MqY^^ f (1 + /i)"^/^esinTt 



Atin-Atout ^ 4sec — 



Vyr; \R^J\RqJ\mJ \{l~b^){l + ecosTt){l-e^f^ 

(7) 

where Tt is the true anomaly at the time of mid-transit and we have neglected the second 
order effect of change in the impact parameter. For a Jupiter-sized planet with a 4 day 
orbit around a solar-sized star, the difference in ingress and egress durations increases to 
3.4 minutes. While Kepler and/or high-precision follow-up observations might be able to 
measure such effects for a relatively small number of bright stars, we expect other methods 
will typically provide more powerful constraints on the eccentricity. 



2.3. Uncertainties in Stellar Mass and Radius 



For many transiting planets photometric observations will be so precise that uncertain- 
ties in the stellar properties will limit the accuracy of fo or ff, and hence the eccentricity 
constraints (unless there are significant constraints from RV observations or time of sec- 
ondary transit; see §3.41) . From the perspective of stellar modeling, the physical properties 
of a star (e.g., radius, and density) are a function of at least three key variables: stellar 
mass, composition, and age. Hence stellar modeling can only provide powerful constraints 
on the star's physical properties if there are at least three observational constraints. The 
effective temperature (Tes) and metallicity {[Fe/H]) can be accurately derived from a sin- 
gle high-precision spectroscopic observation. When a precise parallax (vr) is available, the 
stellar luminosity (L^) can be calculated from the apparent magnitude (V; and a bolometric 
correction). The three constraints {Tes, [Fe/H], L^) can be combined with stellar evolu- 
tion tracks to determine the stellar mass, radius, density, etc. (e.g.. Ford et al. 1999), with 
some well-known degeneracies (e.g., hook region). For relatively bright and nearby stars 
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in the California and Carnegie Planet Search and Hipparcos catalog, this method has been 
used to estimate stellar parameters such as the mass and radius (Valenti & Fischer 2005; 
Takeda et al. 2007). We applied the Bayesian stellar parameter estimation code of Takeda et 
al. (2007) to calculate density directly from the joint posterior distribution, accounting for 
correlations between parameters and the non-linear transformation between observable and 
derived quantities. We find that this method can determine the mean stellar density with 
an accuracy of ~ 4 — 10%, with random uncertainties dominated by the parallax. 

We caution that there may also be systematic uncertainties due to the stellar cvohi- 
tionary tracks. Such systematic effects could lead to or being systematically over- or 
under-estimated and that the error could be highly correlated with the stellar type or other 
stellar parameters. Ideally, such systematics would be mitigated if the densities of stellar 
models were validated by independent observations such as double-lined eclipsing binaries for 
similar type stars. Both the CoRoT and Kepler missions will contribute to stellar astronomy 
and arc likely to contribute towards improving and testing the precision of stellar models. 
In the absense of externally validated models, systematic uncertainties in stellar models may 
limit this techinque for systems with small eccentricities. Fortunately, the stellar density en- 
ters only to the one third power, so potential systematic effects are unlikely to be significant 
for most planets that have sizable eccentricities. In practice, the eccentricity analysis should 
be coupled to a sensitivity analysis that tests whether conclusions are sensitive to the choice 
of stellar models. 

Unfortunately, many target stars for Kepler will be too faint and distant to have well 
determined parallaxes. One alternative approach is to replace with a spectroscopically 
determined stellar surface gravity (logs'). Unfortunately, log (7 is quite difficult to measure 
precisely. For a high-quality spectroscopic observation (e.g., ax^^ ~ lOOK and crpe/H] ^ 
O.ldex), the formal uncertainty in log 51 would typically be ~ O.ldex. Unfortunately, there 
are typically significant correlations with both T^s and [Fe/H] (Valenti & Fischer 2005). 
For a main sequence solar-type star, the uncertainties in the atmophseric parameters would 
translate to an uncertainty in pj (and hence r^) of ^ 8%. In some cases, additional 
constraints on stellar properties such as rotation rate, strength of Ca HK emission, presence of 
Li, and/or asteroseismology may be able to further constrain stellar parameters and improve 
the stellar sensitivity for measuring eccentricities. 

3. Statistical Methodology 

In this section, we outline a few of the statistical approaches that could be applied 
to translate measurements of Tq and/or Tf, into eccentricity constraints. First, we outline 



- 10 - 



a Bayesian approach to calculating the joint posterior probability distribution for the ec- 
centricity and argument of periastron for individual planets. These constraints would be 
particularly valuable when combined with additional constraints on dynamical properties of 
the system. Second, we show that the distribution of normalized transit durations contains 
significant information about the eccentricity distribution of a population of planets. 

3.1. Individual Planets 

For planets with high signal-to noise transit light curves and well-measured stellar prop- 
erties, we suggest a Bayesian framework for determining the constraints on the eccentricity 
and argument of pericenter. To illustrate this method, we adopt non-informative prior proba- 
bility distributions of p(e) = 1 for < e < 1, p{lj) = l/(27r) for < a; < 1, andp(cosi) = 1/2 
for — 1 < cosi < 1. We assume that the observed u is normally distributed about the true 
value of Tb (i.e., ff, ~ N{ti,, cr^J) with a^^ reflecting the uncertainty in the stellar parameters. 
In Fig. [HI we show examples of the joint posterior probability distribution for six possible 
values of ft, . Note that in this figure, u is measured relative to the line of sight, not relative 
to the plane of the sky, as is typical when using radial velocity observations. For actual 
systems, Markov chain Monte Carlo based simulation methods (Ford 2005, 2006) could be 
used to calculate posterior probability distributions allowing for parameters with correlated 
uncertainties and/or non-trivial error distributions (e.g., Holman et al. 2006; Burke et al. 
2007; Takeda et al. 2007), as well as eccentric orbits. 

3.2. Characterizing the Eccentricity Distribution of Transiting Planets 

We have performed Monte Carlo simulations to calculate the expected distribution of 
transit durations for various eccentricity distributions. Here, we consider two limiting cases 
corresponding to low and high signal-to-noise observations: a) only the transit depth, dura- 
tion, and orbital period are measured, so that the observed light curve provides no informa- 
tion about the impact parameter (Fig. [7j), and b) the time of all four points of contact are 
measured so that the observed light curve provides a good estimate of the impact parameter 

(Fig. ED. 
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3.2.1. Impact Parameter Not Measured 

While it would be most desirable to measure the eccentricity for each individual planet, 
this will not always be possible. For some planets (particularly those around faint stars), the 
light curve will not be measured with sufficient precision to measure the orbital inclination 
and sufficient high-precision RV observations will be impractical. In these cases, it will still 
be possible to constrain the mean eccentricity of a population of such planets. Since these 
challenges will be most common for faint stars, we expect that missions such as CoRoT and 
Kepler will detect a large number of such planets available for statistical analyses. 

We consider the distribution of Tq = tD/tD,o, where t/j is the actual transit duration, 
tD^o is the transit duration expected for the same planet and star, but assuming a circular 
orbit (e = 0) and central transit (6 = 0). The distribution for the impact parameter, b, is 
determined by assuming that the inclination is distributed isotropically (subject to the con- 
straint that a transit occurs). We also assume that the observed fo is normally distributed 
about the true value of Tq (i.e., fo ~ N{to, cr'^J) with o"t-„ — 0.05 — 0.15 due to uncertainty in 
the stellar parameters (depending on the stellar properties and available follow-up observa- 
tions). In Fig. [71 we show the distribution of To for six eccentricity distributions. In Fig. [9] 
(right), we plot the mean, standard deviation, and skewness of Tq for an ensemble of planets 
with a given eccentricity. For planets with a circular orbit and an isotropic distribution of 
inclinations. To will have a mean of ~ 0.79 and a standard deviation of ~ 0.22 (Fig. [71 panel 
a). If we assume a uniform distribution of eccentricities (e ~ f/[0, 1)), then the mean Tq 
decreases to 0.64 and the standard deviation increases to 5.2. Alternatively, if we instead 
assume an eccentricity distribution similar to that observed for giant planets (e ~ -R(0.3), a 
Rayleigh distribution with Rayleigh parameter 0.3; Juric & Tremaine 2007), then the mean 
decreases to ~ 0.74 and the standard deviation increases to ~ 0.29. Therefore, it would be 
possible to distinguish between all three models with a plausible sample of planets, assuming 
normally distributed errors and asymptotic scalings for the mean of the distribution, even 
without measuring any impact parameters or RVs. If we were to use only the mean value 
of To, then we would be ignoring the variances and shape of the distribution. Instead, a 
Komogorov-Smirnov test can be used to compare the observed distribution of To to the To 
distribution predicted by a given model. 

To determine how large an observed sample is needed to obtain statistically significant 
results, we perform Monte Carlo simulations. For each trial, we generate two samples of 
transiting planets. The first sample represents a hypothetical observed sample of Npi planets. 
The second sample is much larger and is used to calculate the predicted distribution of To 
for a given theoretical model. For both samples, we assume an isotropic distribution of 
viewing angles, i.e., cost ~ U[—l,l) and u ~ ?7[0,27r). For each trial, we test whether a 
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Kolmogorov-Smirnov test can reject the null hypothesis that the two distributions of Tq are 
drawn from a common distribution. For each Npi we perform several thousand trials; we 
increase A'j,; until we find that the null hypothesis is rejected at the 95% confidence level 
for at least half of the trials. We present our results in Table [U Our simulations show 
that a population of planets on circular orbits (e = 0) could be distinguished from either 
e ~ t/[0, 1) (eccentricities distributed uniformly between zero and unity) or e ~ -R(0.3) 
(a Rayleigh distribution with Rayleigh parameter of 0.3) by using measurements of Tq for 
Npi ^ 30 planets, even without measuring any impact parameters or RVs. Distinguishing 
an observed population with eccentricities distributed as e ~ f/[0, 1) from a model with 
e ~ -R(0.3), would require Npi > 117 planets. Distinguishing between an observed population 
with eccentricities distributed as e ~ i?(0.3) rather than a model with e ~ f/[0,l), would 
require Npi ^ 174 planets. The difference is due to the fact that we use a two-sample 
Kolmogorov-Smirnov test and the sample size for the theoretical model is at least an order 
of magnitude greater than the sample size for the sample of "observed" planets. 



3.2.2. Impact Parameter Measured 

When high quality photometry measures the impact parameter {b), then we can analyze 
the distribution of r^, = tD/tofi, where is the transit duration expected for the same 
planet with the same impact parameter, but assuming a circular orbit. If all transiting 
planets had a single eccentricity, then the expected mean value of Tb is given by (r?,) ~ 
a/1 — (Tingley & Sackett 2005) and the expected standard deviation is (cTt-J — (1 — 
e^)^^^\/l — (1 — e^)-*-/^. We show the distribution of Tf, for several eccentricity distributions 
in Fig. [8] and plot the mean, standard deviation, and skewness for sample distributions of tj, 
in Fig. [9] (left). Note that for small eccentricities, the distribution of is strongly peaked, 
making it possible to perform a significant test of the null hypothesis that terrestrial planets 
have nearly circular orbits with a small sample size ( §3.3p . To account for the uncertainty 
in the stellar parameters, we assume cr^^ ~ 0.05 — 0.15 and normally distributed errors and 
asymptotic scalings. Then, it would be possible to distinguish between various models for the 
eccentricity distribution of observed planets with a few tens of transiting planets (see Table 
|2]). We find that the distribution of Tb provides significant constraints on both the mean and 
the width of the eccentricity distribution. For example, Monte Carlo simulations using the 
Kolmogorov-Smirnov test (similar to those described in §3.2.11 but replacing Tq with Tb) show 
that distinguishing between two Gaussian eccentricity distributions with mean eccentricities 
of 0.18 and 0.36 would require Npi > 30 planets. Similarly, distinguishing between two 
Gaussian eccentricity distributions with a common mean of 0.38 and standard deviations of 
0.10 and 0.24 would require Npi ^ 65 planets based on a Kolmogorov-Smirnov test and a 
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95% confidence level. However, it will be very difficult to measure the higher order moments 
of the eccentricity distribution (see Fig. fTOj) . In Table [2], we list the number of planet required 
for several additional pairs of eccentricity distributions. 

3.3. Testing the Null Hypothesis of Circular Orbits for Individual Planets 

We are particularly interested in addressing the question, "How frequently are terrestrial 
planets on nearly circular orbits?" Therefore, we suggest two statistical tests of the null- 
hypothesis that each planet is on a circular orbit. 

3.3.1. Low S/N Light Curves & Long Duration Transits 

For an eccentric orbit with the transit near apocenter, the transit duration can be 
significantly greater than would be expected for the observed orbital period. It is possible 
to measure a minimum eccentricity for long transits, even for low S/N transits and without 
measuring RV parameters. Since a non-zero impact parameter can only decrease the duration 
of transit (relative to a central transit), it will be possible to reject the null hypothesis of 
a circular orbit (and hence detect a non-zero eccentricity) for planets with long duration 
transits such that Tq or Tfe greater than unity. Unfortunately, the geometric probability of a 
transit occurring is greater for short transits when the planet is closer to the star. Therefore, 
if eccentric planets are common, then there will be more transits with r < 1. Fortunately, 
even a small number of long transits can provide significant constraints for the eccentricity 
distribution. 

3.3.2. High S/N Light Curves & Stellar Models 

For short-period planets that can be assumed to have tidally circularized, Tf, = 1, so high 
precision light curve observations may provide additional constraints on the star's properties. 

-1 /o 

For example, Holman et al. (2007) measure pj to ~ 1.5% accuracy for TrES-2 (V=11.4). 
When combined with theoretical models, this allows for very accurate determinations of the 
stellar and planet radii (Sozzetti et al. 2007). This will enable tests of the null-hypothesis 
that a planet is on a nearly circular orbit. If it is accurate, then the stellar properties would 
be measured quite precisely and all measurements should be self-consistent. On the other 
hand, if the planet is actually eccentric, then assuming a circular orbit could result in an 
inconsistency. For example, by comparing the allowed stellar models to a spectroscopic 
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measurement of log (7 one could recognize planets with eccentricities exceeding ^ 0.15 for 
solar-like stars. Alternatively, the putative location of the star in the (Tgfj, p^, [Fe/H]) 
parameter space could be inconsistent with any stellar model (e.g., Sozzetti et al. 2007). 
The sensitivity of this test will vary significantly with Tgg and be most powerful for stars 
slightly cooler than the Sun. Unfortunately, this method is unlikely to be effective for planets 
with orbital periods larger than a few days, since they may have eccentric orbits. 

3.4. Role of Radial Velocity Observations 

The realization that many giant planets are on eccentric orbits was the result of RV 
surveys (e.g., Butler et al. 2006 and references therein). Clearly, RV observations can measure 
orbital eccentricities, provided that there is a sufficient number, phase coverage, and time 
span of observations with sufficiently high precision (Ford 2005). In fact, the observational 
constraints from radial velocity observations and the photometric method that we describe 
are complimentary. Our photometric method is most sensitive for measuring eccentricities 
of planets if the pericenter direction is pointing towards or away from the observer (Fig. [6]). 
On the other hand, the radial velocity method is most sensitive to measuring eccentricities 
of planets with pericenter nearly in the plane of the sky (Ford 2005; Laughlin et al. 2005). 

Unfortunately, many of the CoRot and Kepler target stars (V~ 9 — 16) will be much 
fainter than the typical targets of RV surveys (V~ 5 — 9). Ground-based transit surveys 
have obtained follow-up RV observations for candidate planets, but typically at a relatively 
modest precision. While there are plans to obtain RV confirmation of planet candidates 
identified by CoRoT or Kepler, these observations will require a large investment of tele- 
scope time, especially when attempting to detect Earth-mass or even Neptune-mass planets 
around relatively faint stars. Even when RV observations can confirm terrestrial planet can- 
didates found by CoRoT or Kepler, they will likely make use of the known orbital period and 
phase to observe at near the extrema of the RV curve. Measuring the eccentricity requires 
measuring the shape of the curve and hence many observations spread across a broad range 
of orbital phases. The amplitude of the deviations of the RV curve from that expected for 
the same planet on a circular orbit is less than the total RV amplitude by a factor of e (to 
lowest order in eccentricity). Assuming a single planet with a orbital period known from 
photometry, a detection of eccentricity for a planet with e ~ 0.3 will typically require an 
order of magnitude more observing time than would be required to detect the planet's total 
RV amplitude (at the same level of significance and assuming observations nearly evenly 
distributed in orbital phase). As an additional complication, the RV method measures the 
refiex velocity of the star induced by all planets, not just one transiting planet. If a star 
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harbors a (potentially undetected) non-transiting planet with a RV amplitude comparable 
to or greater than the radial velocity amplitude due to the epicyclic motion of the transiting 
planet, then RV observations would face an even greater challenge in measuring the tran- 
siting planet's eccentricity. Despite these challenges, we certainly encourage radial velocity 
observations for constraining eccentricities whenever they are practical. 

4. Discussion 

We have described how photometric observations can constrain the eccentricities of 
individual transiting planets and characterize the eccentricity distribution of a population 
of planets. For each planet, we can test the null hypothesis that it is on a circular orbit 
and calculate the posterior probability distribution for the eccentricity. A combination of 
such analyses for several transiting planets could be used to characterize the eccentricity 
distribution of a population of planets. For example, this method could be used to investigate 
how the fraction of eccentric orbits varies with the orbital period or physical proprieties of 
the star and planet. We expect this type of analysis to become increasingly powerful given 
the rapidly growing number of known transiting planets. 

This type of analysis will be particularly valuable for low-mass transiting planets or tran- 
siting planets around faint host stars. In both cases, radial velocity follow-up observations 
will be extremely challenging. This will be the case for many transiting planet candidates 
found by space-based transit searches, such as CoRoT and Kepler. The capability of these 
missions to discover terrestrial-mass planets is particularly exciting. Our method could could 
determine if terrestrial planets with low eccentricities like the Earth are common or rare. 
This would provide significant constraints on theories proposed to explain the eccentricities 
of extrasolar planets. For example, Kepler might find many terrestrial planets with low 
eccentricity orbits, suggesting that the mechanisms that excite the eccentricities of giant 
planets are often ineffective for terrestrial mass planets. In this scenario, those terrestrial 
planets that do have large eccentricities might typically be accompanied by nearby giant 
planets, suggesting that it is the giant planets are responsible for exciting the eccentricities 
of terrestrial planets (Veras & Armitage 2005, 2006). Alternatively, Kepler might find that 
terrestrial planets are much more common than giant planets, and yet they still commonly 
have large eccentricities. This could arise from eccentricity excitation mechanisms that do 
not require giant planets, or due to interactions with previous giant planets that have since 
been ejected, accreted or destroyed by the star (e.g.. Ford et al. 2005; Raymond et al. 2006; 
Mandell et al. 2007). 

Our method also provides a means for studying the tidal evolution of short-period 
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planets. Tidal effects are likely to circularize planets with sufficiently short orbital periods. 
For short-period planets, we can compute the tidal circularization timescale (tdrc) based on 
the properties of the star and planet. If there is a sharp transition between circular and 
eccentric orbits, then this could be used to place constraints on tidal theory (e.g., Zahn & 
Bouchet 1989; Melo et al. 2001; Mathieu et al. 2004). For eccentric planets with relatively 
short tcirc, it may be possible to place a lower limit on the Q factor the is related to the 
planet's internal structure (e.g.. Ford et al. 1999; Bodenheimber et al. 2001; Maness et al. 
2007; Mardling 2007). 

In a Bayesian framework, one can calculate the posterior probability distribution for 
the fraction of each orbit during which the planet-star separation is between an inner and 
outer cut-off. If the cut-offs are set to be the putative boundary of the habitable zone, then 
we can then ask, "What fraction of terrestrial planets are in the habitable zone for some/at 
least half /all of their orbit?". The results of such investigations could have implications for 
the climates of potentially habitable planets, the frequency of such planets, and the design 
of future missions that aim to detect and characterize nearby Earth-like planets that could 
harbor life (e.g., Marcy et al. 2005). 

We have demonstrated that it is practical to collect sufficient photons to characterize the 
eccentricity distribution of terrestrial extrasolar planets, but we assumed that limb- darkening 
parameters can be well constrained by some combination of stellar modeling and external 
observations. Our analytical estimates have not incorporated limb darkening effects or poten- 
tial systematic uncertainties in stellar models. Future search should address both of these 
effects. Multi-wavelength observations (particularly in the infrared) could be particularly 
useful for addressing both these issues. In particular, we plan to investigate the potential for 
combinations of space-based detections and ground-based follow-up observations to improve 
the characterization of the eccentricities of transiting planets. 

Finally, we note that this method for characterizing the eccentricities of terrestrial plan- 
ets from transit light curves underscores the importance of developing and validating precise 
and accurate stellar models. Uncertainties in stellar parameters models are expected to 
dominate the error budget for bright target stars. The potential for systematic uncertainties 
due to stellar modeling will make it particularly challenging to study low eccentricity sys- 
tems as a function of stellar properties. Fortunately, we these concerns would not preclude 
our method from being applied to terrestrial-sized planets recognizing the relatively large 
eccentricities typical for giant planets. 

We thank G. Bakos, D. Charboneau, M. Holman, D. Latham, G. Takcda, and an anony- 
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Table 1. Number of Planets Required to Distinguish between Distributions (Low S/N) 
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^We list the number of planets required to distinguish an observed distribution (taken from the top row) from a 
theoretical distribution (taken from the left column) at the 95% confidence level based on a Kolmogorov-Smirnov 
test applied to Tq. 



Table 2. Number of Planets Required to Distinguish between Distributions (High S/N) 
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'^Wc list the number of planets required to distinguish an observed distribution (taken from the top row) from a 
theoretical distribution (taken from the left column) at the 95% confidence level based on a Kolmogorov-Smirnov 
test applied to fft with =0.15 (bottom) and 0.05 (top). 
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Fig. 1. — Geometry of Transit: Here we illustrate the path (dashed hne) of a planet (small 
circles) as it transits a star (large circle) . The vertical dotted lines connect the planet to the 
schematic light curve (solid curve) at the time of each of the four points of contact. Measuring 
the times of all four points of contact allows for a measurement of the impact parameter 
(6) and significantly increases the precision of the constraint on the planet's eccentricity. If 
only the total transit duration {tu) is measured, then the eccentricity distribution can still 
be constrained for a larger population of transiting planets. 
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Fig. 2. — Transit Duration as a function of eccentricity (e) and argument of periastron (uj): 
The vertical axis shows r^,, the ratio of the actual transit duration to the transit duration 
for a similar planet on a circular orbit (with the same sizes, orbital period and impact 
parameter). The curves show the analytic approximation for the transit duration (Eqn. 1) 
for four values of eccentricity, 0.25, 0.5, 0.75, 0.9. The points are the exact durations for a 
planet on a Keplerian orbit. For a planet on a circular orbit, — 1, but for eccentric planets 
Tb can be larger (for planets that transit near apocenter) or smaller (for planets that transit 
near pericenter). Thus, a measurement of the transit duration can be used to constrain a 
combination of e and cu (measured from the direction of the observer), as well as to place a 
lower limit on the eccentricity. 
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Fig. 3. — Probability Distribution for Transit Duration at a given eccentricity (e): Here we 
show the probabihty distribution (sohd histogram) for ti, for an ensemble of planets with a 
single fixed eccentricity, assuming a uniform distribution of the argument of pericenter and 
an isotropic distribution of inclinations. The dotted curves are the cumulative probability 
distributions. For large eccentricities, there is a small fraction of very long transits (t;, > 2 
off the scale) . 
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Fig. 4. — Transit light curves illustrating effects of limb-darkening. Here we show the 
observed intensity {I{z); normalized to total stellar flux, 7^) as a function of z, the projected 
separation from the center of the star (measured in stellar radii). The vertical dotted lines 
indicate the points of contact. For this illustration, we adopt a quadratic limb-darkening 
model (Mandcl & Agol 2002) and limb darkening coefficients of 71 = 0.4382 and 72 = 0.2924 
based on ATLAS models for solar-like star in V band (Claret 2000). Top: For a Jupiter-sized 
planet. Bottom: For an Earth-sized planet. 



-27- 



1.2 - 



1 - 




0.8 



1.2 - 



1 - 



0.8 



1 — I — I — I — I — I — ry— I — r 



HHH 1 I I I 



i 



J I I I I I I I L 



I I I I I I I I I I I I I I I I I 



H- 



++++ 



++++ 



I I I I I I 



I I I 1 1 1 1 



I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 







++++ 



++++ 



++++■ 



-A ?o 

a 

II 

ro 



1 1 1 1 1 1 1 1 1 1 1 



1 1 1 1 



0.5 
b 



0.5 1 -0.5 0.5 



U1 + U2 



Fig. 5. — Correlations between transit duration and impact parameter or limb darkening 
coefficients. Here we shows contour of constant Ax^ = 1, 4, and 9 with in/tD (the ratio of the 
model transit duration to the actual transit duration) on the y-axis. The various coUumns 

have ,T-axcs of the impact parameter {b; left), the sum of the quadratic limb darkening 
coefficients {ui + U2; center), and the difference of the quadratic limb darkening coefficients 
{ui — U2\ right). The top row shows a 1 planet and the bottom row shows a 2R^ planet, 
both assumed to be at 1 AU from a solar mass star. 



- 28 - 




Fig. 6. — Posterior Joint Probability Distribution for e and uj\ Here we consider the eccen- 
tricity constraint based on the measured value of for an individual planet with a measured 
impact parameter. We show contours equivalent to 1, 2, 3, ...-a bounds on the combination 
of eccentricity and argument of pericenter. The shaded regions are excluded by the mea- 
surement of Tfe. The panels correspond to rj, = 0.25 (top left), 0.5 (top center), 0.75 (top 
right), 1 (bottom left), 1.25 (bottom center), and 1.5 (bottom right). Here we assume that 
the detection probability (given that a transit occurs) is independent of transit duration and 
that the measurement of f;, is normally distributed and has a standard deviation, 0"^-^ = 0.1. 
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Fig. 7. — Distribution of Transit Duration for Various Eccentricity Distributions: We show 
histograms (sohd curves) and cumlative distributions (dashed curves) of To, the ratio of the 
observed transit duration to that expected for the same planet, star, and orbital period, but 
a circular orbit and a central transit. Panel a corresponds to only circular orbits, panels b-f 
correspond to a Rayleigh eccentricity distribution with Rayleigh parameter of 0.1 (b), 0.2 
(c), and 0.3 (d), 0.4 (e), and 0.5 (f). 
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Fig. 8. — Distribution of Transit Duration for Various Eccentricity Distributions: We show 
histograms (sohd curves) and cumlative distributions (dashed curves) of r^, the ratio of the 
observed transit duration to that expected for a similar planet, star, impact parameter, and 
orbital period, but a circular orbit. Panel a corresponds to a uniform eccentricity distribution, 
panels b-f correspond to a Rayleigh eccentricity distribution with Rayleigh parameter of 0.1 
(b), 0.2 (c), and 0.3 (d), 0.4 (e), and 0.5 (f). 
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Fig. 9. — Moments of the Normalized Transit Duration Distribution: Here we consider 
several ensembles of transiting planets, each with a single fixed eccentricity. For each eccen- 
tricity, we calculate the distribution of (left; for transits with measured impact parameters) 
and To (right; for transits without measured impact parameters). Here we plot the mean, 
standard deviation, and skewness for both both tj, and To for each eccentricity. We show 
analytic approximations with curves when available. 
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Fig. 10. — Distribution of Transit Duration for Eccentricity Distributions with Common 
Mean and Variance: Top: Histograms of Tq, the ratio of the observed transit duration to 
that expected for the same planet, star, and orbital period, but a circular orbit and a central 
transit. We show results for two eccentricity distributions: normal (solid) and Rayleigh 
(dashed). The mean and variance of the normal distribution have been chosen to match 
that of a Rayleigh distribution with Rayleigh parameter 0.3. The resulting distributions of 
To are so similar that it would be extremely difficult to distinguish between these eccentricity 
distributions. Bottom: Same as above, but for histograms of Tb, the ratio of the observed 
transit duration to that expected for a similar planet, star, impact parameter, and orbital 
period, but a circular orbit. Again, the distributions of are too similar to distinguish 
between the eccentricity distributions. 



